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Abstract 

This article demonstrates an undergraduate experiment for the determination of specific heat 
capacities of various sohds based on a calorimetric approach, where the soHd vaporizes a measurable 
mass of liquid nitrogen. We demonstrate our technique for the metals copper and aluminum, the 
semi-metal graphite and also present the data in relation with Einstein's model of independent 
harmonic oscillators and the more accurate Debye model based on vibrational modes of a continuous 
crystal. Furthermore, we elucidate an interesting material property, the Verwey transition in 
magnetite occurring around 120-140 K. We also demonstrate that the use of computer based data 
acquisition and subsequent statistical averaging helps reduce measurement uncertainties. 
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I. THEORETICAL MOTIVATION 



The departure of the heat capacity from the classical Dulong-Petit law at low tem- 
peratures is one example of the success of quantum mechanics in describing experimental 
observations.- The idea highlights the confluence between apparently dissimilar phenomena 
such as Planck's formula used in studying blackbody radiation, quantization of atomic vi- 
brations, macroscopic heat capacities and propagation speeds of waves. An experiment that 
is accessible to undergraduates for measuring low temperature heat capacities and correlat- 
ing results with Einstein's and Debye's quantized descriptions, is therefore, of fundamental 
significance. We first quickly refresh some of the underlying theory. 

A. Einstein's model 

The details on the specific heat capacity of solid can be looked up in any standard text 
on solid state physics.^ In summary, Dulong and Petit, using the equipartition theorem of 
classical thermodynamics, showed that the molar heat capacity of metals has a constant 
value of 3R, R being the molar gas constant. However, as we move to lower temperatures, 
the heat capacity shows significant variation. This deviation inspired Einstein do the first 
quantum mechanical calculation of the specific heat capacity.-*^ He assumed a solid compris- 
ing independent, three-dimensional oscillators (per unit mole), all possessing the same 
fundamental frequency / and derived the following temperature (T) dependence of the molar 
heat capacity Cy, 



where ks is Boltzmann's constant and the variable 6e = hf/ks, with dimensions of temper- 
ature, is called the Einstein temperature. In the high temperature approximation, T >> 6e, 
Cy approaches SNks = 3R as predicted by the Dulong-Petit law. 

B. Debye's model 

Peter Debye considered^ the vibrational modes of a continuous medium, as opposed to 
the vibrations of independent atoms predicated in Einstein's model. The density of states 
g{f) for a three dimensional solid is 47r/^/c^, where c is the propagation speed of the wave 




(1) 
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inside the solid. Since the speeds vary for the longitudinal (q) and the doubly- degenerate 
transverse (q) waves, one may also write, 

g(f)df = 47rf(\ + ^)df. (2) 

Using the above frequency spectrum and Planck's formula, the total vibrational energy of 
the crystal is, 

E = ^(i_ + 4)TM -^dx. (3) 

Here we have made the substitution x = hf/ksT. The Debye frequency fu is a cut-off value 
restricting the total number of modes to 3A^, similar to the atomic vibration case. This is 
achieved by letting, 



f%{f)df = 3N. (4) 

^0 



fl-'4(^^^] . (5) 



The cut-off procedure yields. 

Often, one defines the Debye temperature 6d = hfij/ks, which is a useful parameter as it 
determines the density of oscillators (as defined earlier) and the wave propagation speeds. 
The former can be used to predict the mass density of the material. The upper limit in the 
integral. Equation (|3]), therefore becomes Oo/T. Differentiating with respect to temperature 
yields the expression. 



C. Numerical predictions 

Based on either Einstein's model. Equation ([1]), or Debye's model. Equation ([6]), one 
can make numerical predictions of heat capacities. Einstein's model is straightforward to 
simulate, while a numerical integrator is required for Equation (|H]). For example, using 
Mathematica, the integral avaluates to, 

dx = -^-^+yHn{l- e^) + 3y' Lz^le^) - Qy Ltsie^) + 6 Lz4(e^), (7) 
Jo - 1 15 4 

where y = 6d/T and the polylogarithm function^ is defined as Lin{z) = H'^^^z^ /p^. The 
numerically predicted values can be graphed and compared with experimental results. 
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II. THE EXPERIMENT 



A. An overview 

In this work, we present an experiment for the thermal physics laboratory that aims at 
experimentally determining low temperature heat capacities, down to ~ 100 K, (a) showing 
deviation from Dulong- Petit law, (b) a comparison with numerical predictions of Einstein's 
and Debye's models and (c) an estimation of 9e and Od. In addition, the heat capacity is di- 
rectly linked to material properties and provides insight into (d) structural (re) arrangements 
of atomic constituents of materials. For example, the specific heat capacity changes near 
phase transitions,- such as at the Curie point in ferromagnets.- It's quite an enriching ex- 
perience for students to investigate these properties inside the laboratory setting. 

Thompson and White have presented^ details of a beautiful experiment measuring the 
latent heat of vaporization of liquid nitrogen and the specific heat capacity of various metals. 
Their method is based on calorimetric heat exchange between the solid and liquid nitrogen. 
The present discussion is a straightforward extension of their work, the added feature being 
the determination of low temperature heat capacities enabled by a statistical minimization 
of uncertainties — a lucid example of using statistics to one's advantage. We show results 
from experiments performed on Cu, Al, and graphite and go on to illustrate the anomaly in 
the heat capacity in a ferrimagnetic material at the so-called Verwey transition. 

B. Experimental scheme 

The arrangement of the apparatus is shown in Figuredl A force sensor with an accuracy of 
0.01 N (1 g) (Vernier Instruments DFS-BTA) is interfaced to the computer and measures the 
weight (and hence mass) of a styrofoam cup containing the continuously vaporizing nitrogen. 
The solid, whose heat capacity is to be determined, makes contact with cold nitrogen vapor 
from the boiling liquid placed inside a vacuum flask which serves as our vapor cryostat. 
The temperature of the solid is monitored by a silicon diode that is physically secured 
with teflon tape while thermal grease ensures uniform thermal contact between sensor and 
solid. The level of the liquid nitrogen is kept constant by continuous refilling through a 
funnel. As the desired initial temperature of the solid, Ti, is achieved, it is swiftly dropped 
into the styrofoam cup containing liquid nitrogen. A rapid hissing sound and effervescence 

4 



ensues, partly due to the Leidenfrost effect^ and the cup jostles. Finally, the rapid movement 
settles and the background rate of vaporization is re-established. In equilibrium, the solid's 
temperature is 77 K. A representative variation of the solid's temperature during the course 
of the experiment is illustrated in Figure |5J 



C. Measuring temperature 

We use an ordinary off-the-shelf silicon diode as our cryogenic temperature sensor.— 
The underlying principle of the temperature measurement is the diode equation. If = 
Io{T, Eg)(exp {eVf /kBT) — l), where // is the forward-biased current through the diode, 
Vf is the forward-biased voltage, ks is Boltzmann's constant, T is the temperature and /q 
is the reverse saturation current, which itself depends on temperature and the band-gap Eg. 
It can be shown^^ that the equation relating the temperature with diode voltage is given by, 

Vf{T) = ^ - (log(a) + I log(T) - log(/,))^, (8) 

with, 

(. 3 

and q, m, A and r are, respectively, the charge and mass of the electron, cross- sectional 
area of the diode junction and the momentum scattering time.^ A plot of Vf versus T is 
approximately linear in our temperature range of interest, i.e., from ~ 100 K to room 
temperature. 

To provide a constant forward biased current // = 10 fiA, a current source can be easily 
built using an operational amplifier (TL081), Zener diode (2.7 V), and resistors. The circuit 
is shown in Figure El The resulting voltage Vf is directly read into the computer fitted 
with a data acquisition card (National Instruments PCI-6221) and a Labview programme 
automates the data acquisition. 



D. Why not to use a thermocouple? 

Note that unlike the wire thermocouple, the diode gives a stabler and more accurate 
measurement for the sample held inside the nitrogen vapor. The thermocouple operation 
is principled after the Seebeck effect, wherein a temperature gradient along the length of a 
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conductor results in an emf. The thermocouple wire protruding above the liquid nitrogen 
surface in the cryostat is in fact placed inside a spatially extended temperature gradient. 
Therefore, the induced Seebeck voltage originates from the entire thermocouple wire and 
not only the welded tip making contact with the solid surface. Hence thermocouple mea- 
surements don't represent the temperature of the solid. 

E. Surface or bulk temperature? 

An additional concern that may arise is whether the surface temperature truly represents 
the bulk temperature. Typically our solids are cylinders of diameter 12 mm and length 
32 mm. The Si diode indeed measures the surface temperature of the solid, but we can 
make an intelligent guess about the internal temperature based on the Biot number, Bi = 
HLf./ K, where H is the coefficient of heat transfer in the presence of nitrogen vapor, is 
the ratio of the solid volume to its surface area and k is the thermal conductivity. The value 
of H depends on the temperature of the nitrogen about its boiling point,— but for a gap 
greater than 40 K, it is approximately 3000 W/(m^ K). A simple calculation shows that for 
the Al and Cu samples, the Biot numbers are approximately 0.03 and 0.016, both of them 
are smaller than one, indicating that the temperature is uniformly distributed within the 
solid. 

F. Data analysis and statistical minimization of errors 

Data from a representative experiment are shown in Figure HI For the data processing 
part, we start off by constructing straight lines as best fits, let's call them a and b to the 
background evaporation rates before and after immersing the solid. The vertical displace- 
ment between these lines is a measure of the additional nitrogen vaporized by the heat 
fiowing from the hotter solid. Am. The specific heat capacity is determined from the change 
in mass. 



where Umoies is the number of moles of the solid and is the latent heat of vaporization of 
nitrogen. For the temperature dependence of C^, one has to simply repeat the experiment 
by varying the temperature Ti and calculating the corresponding decrease in mass Am. 




J/(mol K), 



(10) 



6 



Figure [5] shows the experimental results from solid copper, acquired at the initial tem- 
peratures Ti of 200 K (subfigures (a) and (c)) and 120 K (subfigures (b) and (d)). The 
vertical lines are the error bars, let's call them Um, arising from the instrumental least 
count, = 1 g- For a meaningful analysis, we strictly require Am > Um- In Figure a), 
this condition is easily satisfied. The lines prior and after immersion, when extrapolated, 
do not statistically overlap. The situation, however, changes when the initial temperature 
of the solid Ti decreases, resulting in smaller values of Am. This trend is sampled in Fig- 
ure [5t^b): the instrumental uncertainties overwhelm the change in mass making it impossible 
to obtain a statistically reliable value of Am. 

The statistical work-around this limitation lies in appreciating that we are using the 
entire lines prior and after the immersion of the solid for the determination of Am, and 
not just the terminal values at the juncture of the transferral of the solid. Therefore, we 
are interested in the errors in the lines, vis-a-vis, the uncertainties in the slope Ug and the 
intercept Ui. These ensemble uncertainties can be computedr^ using, 



where di is the deviation of the i'th experimental point from the corresponding point on the 
best-fit line, m is the average mass, D is the sum of squares of deviations and n is the total 
number of points in each line. Computer-based acquisition generates large amounts of data, 
increasing n and hence reducing the uncertainties in the measurement. Based on Ug and Ui, 
we can draw confidence bands, with larger ra's resulting in even tighter bands. These bands 
are illustrated in Figure [S](c) and (d) by a trio of lines: the middle line is the curve of best 
fit and the top and bottom lines represent the extrema of the band. 

The best estimate for the reduction in mass is Am and its uncertainty UAm is based on 
the maximum (Ami) and minimum (Am2) differences. These maximum and minimum are 
defined through the inset of Figure [5](e) . The uncertainty is. 



This work-around ensures that tiAm < Am even in cases where Um > Am. A numerical 
example is in order here. For an initial solid temperature of 120 K, the individual mass 
uncertainty Um = 1 g is larger than the change in mass Am = 0.72 g, but the uncertainty 




(11) 



(12) 




(13) 
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calculated after the statistical averaging procedure is UAm = 0.43 g, which is smaller than 
Am. Finally, the uncertainty in C^, uc^ is inferred from using the well known error 
propagation formulas,— 

uc^ = — — , (14) 

nmoiesy-^ 1 — I I ) 

showing that smaller initial temperatures yield higher uncertainties mc„. This trend is di- 
rectly observable in our results. 



III. RESULTS AND DISCUSSION 
A. Einstein fits for Cu and Al 



Results for the metals Cu and Al are shown in Figure [61 The data is numerically fitted 
to the Einstein curve. Equation ([T]), yielding 6e^s of 278 and 284 K, respectively, agreeing 
reasonably well with the nominal values {e.g., 248 and 306 K).— Students can also observe 
that the heat capacity for Al is smaller than Cu, implying higher Einstein temperatures and 
a higher frequency /. Students can get an appreciation of this fact by recognizing that the 
Al atom is lighter than Cu (smaller m) and moreover, the elastic constant is higher (higher 
k) (volume expansivity is lower). This means that for Al, the ionic frequency (/ oc ^Jkjm) 
will be higher, resulting in a higher Einstein temperature Qe = hf/ks. It is clear that the 
data obtained from this simple and inexpensive route show considerable agreement with 
qualitative predictions. 



B. Debye fits for Cu and Al 

The experimental data for Cu and the Debye fit with 6^ = 350 K is shown in Figure [Tl^a), 
with the agreement being exceptionally good at lower temperatures. The best estimate of 
On is found by minimizing the variance between the experimental and numerical values, 
(C^P^-d°""^)',the process is illustrated in Figure [Tl^b). Our result is in excellent 
agreement with the published value of 343 K.^ Analogous results from Al yield a Debye 
temperature of 450 K, compared to the literature value of 428 K.- 
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C. Heat capacity of graphite and Verwey transition in magnetite 

The specific heat capacity for graphite is shown in Figure [8](a) showing good agreement 
with pubhshed work.— A commonly available ferrite is magnetite, Fe304, well known for its 
commercial usage in ferrite cores of inductors and transformers, and its interesting magnetic, 
electrical and structural properties. For example, magnetite has the iconic inverse spinel 
structure, important for the description of high temperature superconductivity in copper 
oxide based ceramics. For the present purpose, we are interested in highlighting an anomaly 
in magnetite's heat capacity — the Verwey transitioi>i^ which occurs around 120-130 K. At 
temperatures above the Verwey transition, the ferrite is metallic and electrically conducting 
while at lower temperatures it is an electrical insulator, its conductivity decreasing by several 
orders of magnitude. Besides, the change in electrical conductivity, the heat capacity also 
shows anomalous behaviour. Our experimental results are shown in Figure Mjo) and are in 
agreement with Parks's original experiment^ in 1926. 

In conclusion, the present article gives a practical illustration to determining the specific 
heat capacity at low temperatures using a vapor cryostat for lowering temperature, a silicon 
diode as a temperature probe and a gravimetric technique to track the vaporization of 
nitrogen with and without a solid specimen. The experimental arrangement, which is simple 
yet elegant, can bring home important concepts in statistical mechanics, materials physics 
and provide a direct demonstration of quantization of atomic or crystal vibrations. The 
present work builds on other foundational experiments or numerical techniques reported in 
the same journal.-i^i^ The temperature dependent measurements are easily extensible to 
investigating complementary material properties such as electrical conductivity and band 
gap,— as well as latent heats of phase transitions (such as the at Curie or Neel transitions 
and the superconducting to normal transition of high temperature cuprate superconductors). 
We hope this experiment will be a useful addition to the thermal physics laboratory. 



* Electronic address: sabieh@lums . edu. p kl 

^ C. Kittel, "Introduction to Solid State Physics", (John Wiley & Sons, Delhi, 2005). 
^ A. Einstein, "Planck's theory of radiation and the theory of specific heat", Ann. Phys. 22, 
180-190 (1907). 



9 



^ E. Lagendijk, "A simple determination of the Einstein temperature", Am. J. Phys. 68 (10), 
961-962, (2000). 

^ P. Debye, "Zur theorie der specifischen warme", Ann. Phys. 39, 789-839, (1912). 

^ D. H. Bailey, P. B. Borwein and S. Plouffe, "On the rapid computation of various polylogarithmic 

functions", Math. Comp. 66, 903-913, (1997). 
^ R. H. Dittman and M. W. Zemansky, "Heat and Thermodynamics", (Tata McGraw Hill, New 

Delhi, 2007). 

^ J. M. D. Coey, "Magnetism and Magnetic Materials", (Cambridge University Press, United 
Kingdom, 2010). 

^ C. W. Tompson and H. W. White, "Latent heat and low-temperature heat capacity for the 

general physics laboratory". Am. J. Phys. 51 (4), 362-364, (1983). 
9 F. L. Curzon, "The Leidenfrost phenomenon". Am. J. Phys. 46, 825-828, (1978). 
^ P. R. N. Childs, J. R. Greenwood, and C. A. Long, "Review of temperature measurement", 

Rev. Sci. Listr. 71, 2959-2978 (2000). 
^ J. W. Precker and M. A. de Silva, "Experimental estimation of the band gap in silicon and 

germanium from the temperature voltage curve of diode thermometers". Am. J. Phys. 70, 

1150-1153 (2002). 

^ A. Khalid and M. S. Anwar, "Superconducting Quantum Interfer- 
ence Devices". This technical monograph can be downaloaded from 
http://physlab.lums.eud.pk/ /index.php/Experiments jn_Lab-n 

^ F.P. Incropera and D.P. Dewitt, "Fundamentals of Heat and Mass Transfer", (Wiley, New 
Jersey, 2001). 

^ T. Jin, J. Hong, H. Zheng, K. Tang, and Z. Can, "Measurement of boiling heat transfer coeffi- 
cient in liquid nitrogen bath by inverse heat conduction method", J. Zhejiang Univ. Science A 
10, 691-696 (2009). 

^ G. L. Squires, "Practical Physics", (Cambridge University Press: New York). 
^ H. Ledbetter, "Thermal expansion and elastic constants", Int. J. Thermophys. 12, 637 (1991). 
^ K. Komatsu and T. Nagamiya, "Theory of the specific heat of graphite", J. Phys. Soc. Jpn. 6, 
438-444 (1951). 

^ B.D. Cuhity and CD. Graham, "Introduction to Magnetic Materials", (IEEE Press and Wiley, 
New Jersey, 2009). 

10 



J. Garcia and G. Subi'as, "The Verwey transition — a new perspective", Condens. Matter J. 
Phys. 16, R145-R178 (2004). 

G. S. Parks and K. K. Kelley, "The heat capacities of some metalhc oxides", J. Phys. Chem. 
30, 47-55 (1926). 

Christropher G. Deacon, John R. de Bruyn, and J. P. Whitehead, "A simple method of deter- 
mining Debye temperatures". Am. J. Phys. 60 (5), 422-425, (1992). 



Nitrogen level 

is maintained Force sensor 
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FIG. 1: The experimental arrangement. The solid is allowed to thermally equilibrate with the cold 
nitrogen vapor inside a cryostat and is then swiftly migrated and dropped into a cup holding liquid 
nitrogen. The weight of the cup is constantly monitored while temperature is measured by passing 
a fixed current through and measuring the voltage drop across a Si diode. 
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FIG. 2: Temperature profile of the solid during the course of the experiment: a: represents the 
point at which the solid achieves the desired temperature inside the vapor cryostat; b: shows the 
instant when the solid is lifted from the cryostat; c: is in the region of the decreasing temperature 
during the solid's transit from the cryostat towards and subsequent immersion into the cup, and 
finally d: is the point at which the solid eventually equilibrates at 77 K. 



-12V 
p 




-12V 



FIG. 3: A home-built circuit for the current source to provide 10 /iA current through the silicon 
diode. The current is adjusted by the setting of the variable resistor. 
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FIG. 4: Mass of the vaporizing nitrogen versus time: a: (blue line) background evaporation rate 
of boiling liquid nitrogen before immersing the solid; b: (red line) reestablished evaporation rate 
after the solid is immersed. The mass of the solid has been subtracted in plotting b. 
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FIG. 5: (a) shows the temporal profile of the evaporating mass for Ti = 200 K and (b) the profile 
for Ti = 120 K. The time-stamped measurements of mass and the individual uncertainties in each 
measurement Um are presented as well as the best line fits prior and after the immersion of the 
solid specimen. Note that for ease of visualization, most of the sample points have been removed, 
(c) and (d) show the computed bands based on the uncertainties in the entire lines. The individual 
points and their uncertainties have been suppressed in these latter subfigures. For Ti = 120 K, 
we have Ami = 5.10 g, Am2 = 4.50 g, Am = 4.80 g, UAm = 0.45 g and uc„ = 0.44 J/(mol K) 
and for Ti = 200 K, we compute, Ami = 1-10 g, Am2 = 0.34 g. Am = 0.72 g, UAm = 0.43 g and 
uc^ = 1.4 J/(mol K). The inset (e) illustrates a close-up view of typical measurement bands in 
the vicinity of migrating the solid from the cryostat to the styrofoam cup. The labeling a-f is used 
for defining Am, Ami and Am2 in the text (see Equation I13p. 
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FIG. 6: Experimentally determined heat capacities of (a) copper and (b) aluminum along with the 
Einstein fits. 
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FIG. 7: Experimentally determined heat capacities of (a) copper along with the best fit according 
to Debye's model. The fitting procedure is illustrated in (b) where the variance between the 
experimental points and the numerical estimates is minimized to find the best estimate oi 9d- 
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FIG. 8: Experimentally determined heat capacities of (a) graphite and (b) magnetite, the latter 
highlighting the change in specific heat capacity in the region of the Verwey transition. In (a), our 
experimentally determined values are shown by the squares (■) and the circles (•) are the data 
points from Komatsu's workJ^ 
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